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Abstract —  Natural  decay  of  pollutants  is  rarely  consider  in 
mass  transport  equation.  This  paper  develops  a 
mathematical  model  of  a  filter  based  on  iron  oxide  coated 
sand  for  the  removal  of  chromium  in  groundwater  and 
surface  water.  The  natural  decay  of  chromium  was 
analyzed  using  a  zero  order  decay  reaction.  The  advection 
dispersion  equation  for  this  purpose  is  solved  analytically 
using  the  Laplace  transform  method  and  numerically  by  the 
Fourth  order  Runge  kutta  method  to  determine  the  spatio- 
temporal  distribution  of  the  pollutant  through  the  filter.  The 
control  parameters  are  the  adsorption  coefficient ,  the  initial 
concentration  and  the  chromium  degradation  coefficient.  It 
has  been  proven  that  the  chromium  degradation  coefficient 
strongly  affect  the  concentration  values  of  pollutants  inside 
the  filter.  The  results  of  this  work  should  be  used  to 
manufacture  low-cost  filters  based  on  iron  oxide  coated 
sand. 

Keywords —  Filter,  Sand,  Iron  oxide,  Chromium, 
Adsorption,  advection-dispersion,  Laplace. 

I.  INTRODUCTION 

Water  is  an  essential  component  in  the  universe  and  plays 
an  important  role  in  the  proper  functioning  of  ecosystems 
[1].  Despite  this,  access  to  drinking  water  is  an  increasingly 
recurrent  challenge  worldwide  and  it  is  estimated  that 
nearly  eight  hundred  and  eighty  (880)  million  people  lack 
safe  drinking  water  at  home  [2,  3].  Water  pollution  results 
mainly  from  poorly  planned  urban  and  rural  development, 
industrial  discharges,  agricultural  discharges,  geological  and 
environmental  activities,  domestic  activities  etc.,  it  is  a 
serious  problem  because  it  affects  our  lives  [4].  In  addition, 
the  World  Health  Organization  (WHO)  estimates  that  80% 
of  the  diseases  that  affect  the  population  are  directly  carried 
by  water:  400  million  people  are  permanently  affected  by 
gastroenteritis,  160  million  from  malaria  and  30  million 
onchocerciasis  [5].  The  industry  uses  heavy  metals  in 
various  fields  because  they  have  a  high  electrical 
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conductivity;  however  these  heavy  metals  have  some 
toxicity  to  humans,  leading  in  particular  toxic  effects  on  the 
nervous  system,  carcinogenic  effects  on  the  blood  and  the 
bone  marrow  and  kidney  disorders. 

Chromium  is  one  of  the  most  widely  used  metals  in 
industries,  since  it  has  enough  advantages  for  tanneries, 
textile,  wood  processing,  agribusiness...  Chromium  VI  is 
the  most  problematic  form  of  chromium  since  in  this  form 
chromium  is  very  toxic  and  very  soluble  in  water.  This 
solubility  gives  it  great  mobility  in  ecosystems.  Today, 
following  a  lack  of  compliance  with  laws  or  accidents,  a 
large  number  of  industrial  sites  are  polluted  by  chromium 
Therefore,  the  development  of  low-cost,  easy-to-maintain 
technologies  to  reduce  this  pollutant  in  drinking  water 
sources  to  an  acceptable  level  for  consumption  is  an 
interesting  research  topic. 

The  methods  implemented  today  to  treat  discharges  polluted 
by  chromium  (VI)  are  generally  aimed  at  recovering  this 
metal,  in  order  to  reuse  it  and  thus  reduce  its  eco- 
toxicological  inpact.  In  parallel,  physicochemical  methods 
attempt  to  eliminate  Cr  (VI)  and  electrochemical  methods 
attempt  to  reduce  Cr  (VI)  to  Cr  (III).  Iron  oxide  coated  sand 
has  proven  its  effectiveness  for  the  elimination  of  various 
heavy  metals  (As,  Cd,  Cr,  Pb,  ...)  at  an  elimination  rates 
close  to  100%.  In  this  study,  we  develop  a  mathematical 
model  to  understand  and  predict  the  adsorption  of 
chromium  VI  in  a  fixed  column  bed  made  witchiron  oxide 
coated  sand.  As  Chromium  undergoing  a  natural  decay,  the 
discussion  will  be  focus  effect  of  degradation  coefficient  X 
on  the  spatio-temporal  variation  of  the  pollutants  inside  the 
filter. 

II.  MATERIAL  AND  METHODS 
2.1  Preparation  of  iron  oxide  coated  sand 

The  iron  oxide  coated  sand  was  prepared  using  a  procedure 
similar  to  that  of  Bailey  et  al.  (1992)  [6].  The  author  used 
washed  and  dried  river  sand  with  a  weight  of  200  g  and  a 
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geometric  size  of  0.49  mm  [7].  It  is  obtained  by  mixing  for 
2  min,  80  ml  of  a  solution  of  ferric  nitrate  2M  (Fe  (NO3) 
39H2O)  to  washed  river  sand  and  dried.  He  then  placed  the 
mixture  in  a  drying  oven  at  110  0  C  for  14  hours.  He  finally 
places  the  mixture  in  distilled  water  until  the  runoff  is 
cleaned  at  105  0  C  and  stored  in  closed  bottles  [7]. 

2 .2  Mathematic  al  des  cription 

The  linear  dispersion  advection -equation  is  used  to  describe 
the  transport  of  pollutants  in  this  fixed  column  bed  (Fig.  1). 
The  aquation  was  derived  to  predict  aqueous  concentrations 
of  contaminants  over  time  and  space  in  the  porous  medium 
It  is  written  as  follows :  [8] 


dt  dt  dx2 


dc 

J — 

dx 


(i) 


Where  p  is  the  density,  0  the  water  content,  D  the 
hydrodynamic  dispersion  coefficient,  V:  average  Darcy 
velocity  and  q  the  equilibrium  adsorbed  concentration 
which  is  described  as  follows [8]: 


q  —  kcn  (2) 

Where  k  is  the  adsorption  capacity  and  n  the  adsorption 
intensity. 

In  the  absence  of  the  iron  oxide  coated  sand  which  is  the 
adsorbent,Cr  (VI)  ion  disappeared  naturally  from  the 
solution.  Given  the  fact  that  chromium  is  a  naturally 
decreasing  pollutant,  this  disappearance  rate  was  assumed 
to  be  at  zero  order  as  follows : 


—  del  dt  =  A;  c  =  c0  at  t  =  0  (3) 

Where  X  is  the  zero  order  rate  constant  or  chromium 
degradation  coefficient,  cois  the  initial  Cr(VI)  concentration 
and  c  the  Cr(VI)  concentration  at  the  time  t. 

By  adding  a  decay  term  Ado  the  equation  above  according 
to  Runkel  and  by  considering  a  linear  adsorption  isotherm 
(n  =  1),  after  a  few  transformations,  we  obtain: 


dc  ,  d2c  dc  ~ 

—  =  d  —  -  u  —  -  Ac  (4) 

d  t  dx2  dx 

The  appropriate  boundary  conditions  for  the  model  are 
given  by: 

c  (x  ,0)  =  0  ,  t  >  0; 


c  (0,  t)  =  co,  t  >  0; 

c  (00,  0)  =  0  ,  x  >  0 


(5) 


2.3 Analytical  solution  of  the  mathematical  model 

The  analytical  solution  is  obtained  using  Laplace  transform 
It  is  use  to  transform  one  or  more  partial  derivatives  of  the 
differential  equation  into  algebraic  expressions.  We  call 
Laplace  transform  of  the  function  G  denoted  L  (G),  the 
function  of  complex  variable  P  defined  by: 

L  [G(t)]=  G(p)=J0  e~pt  G(x,  t)dt  (6) 

Where  L  denotes  the  Laplace  transform,  G(t)  is  a  function 
of  time,  G(P) is  a  function  of  the  transformed  space,  it 
corresponds  to  the  Laplace  value.  The  operational  property 


of  the  transformation  that  eliminates  a  partial  derivative  of 
the  equation  is  as  follows : 

L^  =  pG-G(  0)  (7) 

Where  G  (0)  is  the  value  of  G  at  time  t  =  0 
To  reduce  equation  (4)  to  a  more  familiar  form,  we  take: 

u2t 


c(x,  t) 


r(x,  t)exp  -  A'd 

r  2d  4  d 


(8) 


T  N  UX  u2t  -  , 

Let0(x,  t)  = - A  t 

K  J  2d  4  d 


By  substituting  equation  (8)  in  equation  (4),  we  obtain: 

V  =  d^  (9) 

dt  dx2 

The  initial  and  boundary  conditions  (6)  become: 
c  (x  ,0)  =  0  ,  x  >  0 

c  (0,  t)  =Co  exp  (—  +  /Tt),  >  0  (10) 

4  d 

c( oo,  0)  =  0  ,  x  >  0 

The  Laplace  transform  of  equation  (9)  is: 

L(— )  =  d^7  (11) 

KdtJ  dx 2  V  7 

Therefore,  it's  reduced  to  an  ordinary  differential  equation 

The  solution  of  the  problem  is : 

r(x,t )  =  y exp  +  A't]  [ e2£erfc(a  +  ”)  + 

e~2£erfc(a ]  (13) 

Substituting  (13)  in  equation  (8),  the  solution  becomes: 
exp  [(“)  e~2£  erfc  (a  -  +  e2£erfc  (a  +  ^)]  (14) 


Where, a  =  -*=  ,e  =  I (— +  A')  =  ^andF  = 

2 Vdt  ^v4 d  y  l2 VdJ  u2 

Vl  +  2  H 
Finaly: 

j=i{e,P[=a-rt]er/c(^)  +  ^[=a  + 

"5) 

2.4  Numerical  solution  of  the  mathematical  model 

The  numerical  solution  allow  to  take  into  account  more 
comp  lex  heterogeneous  systems,  and  they  can  be  thus  more 
realistic  and  flexible  to  simulate  the  real  conditions.  It  is  not 
very  obvious  to  obtain  an  analytical  solution  of  linear  PDE 
such  as  equation  (4).  However,  there  is  an  approximate 
solution  for  linear  PDEs  obtained  from  different  numerical 
methods.  Therefore,  to  solve  the  EDP  of  pollutant  transport 
in  porous  media,  we  use  numerical  methods.  The  partial 
discretization  method,  where  the  finite  difference 
approximation  is  used  to  obtain  the  discrete  form  of  the 
spatial  derivative.  Then,  the  finite  space  difference  diagram 
in  space  is  used  to  have  discrete  shape  spatial  derivation 
equations  of  equation  (4).  Finally,  the  fourth  order  Runge 
Kutta  system  (RK4)  is  used  to  solve  equation  (4).  This 
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numerical  technique  is  used  to  solve  the  ordinary 
differential  equation  expressed  as  follows: 

%=g(t,C)  (15) 

The  RK4  algorithm  is  written  as  follows : 

ki=At  g(t,  C(t)) 

k2=At  g(t  +^At,  C(t)  4ki) 

k3=  g(t^At,  C(t)  fk2)  (16) 

k4=At  g(t+  At,  C(t)  +k3)  (17) 

C(t  +  At)  =  C(t)H — (ki+2  k2+2  k3+  k4) 

6 

Where,  At  is  the  time  step,  C  (t)  is  the  concentration  of  the 
pollutant  at  time  t,  C  (t  +  At)  is  the  concentration  at  time  (t  + 
At).  The  partial  discretization  method  of  the  finite 
difference  approximation  of  first-order  spatial  derivatives 
is: 


dc 

dx 


i  i 

C-C  . 


Ax 


(17) 


The  discrete  form  of  derivative  of  the  second  order  is: 


d2c 

dx2 


C  -2C.+C 


(18) 


Or 


—  (C]  -C]  )-A'C] 

Ax  V  1  1 


The  indices  (i)  and  (j)  represent  the  discretization  nodes 
along  (x)  and  (t)  respectively.  Ax  is  the  step  of  space.  Thus, 
equation  (4)  can  be  written  in  the  following  reduced  form: 

dC-}  /  j  j  j  \ 

— -  =  g(  t,x,  C  +i,C  ,C-i ) 

Ot  v  iir/ 

,  g{t,  x,  ClX.C'J  =  j~2  (c;+i  -  2C\  +  C'J  - 

(19) 

g(t,x,c;+i,c;,c;j 

=  d^c’+i  -  2C\  +  C-i)  -  d2{c[  -  c'-i) 
-  xc\ 

Withdi.  =—2,d2=  — 

Ax2  Ax 

We  definex  =  iAx ,  t  —  ydtand  (i=0,  1,2...,  Nx)  et  (j=0, 
1,2,...,  Nt) 

The  discrete  version  of  the  initial  condition  associated  with 
equations  (4)  and  (5)  is  expressed  as  follows: 

C°  =  C0  if  i  =  0  (20) 

C°  =  0  if  i  *  0 

Equations  (20)  express  a  column  without  initial 
contamination.  The  discrete  form  of  the  boundary 
conditions  is:C.°  =  C0,C0is  the  initial  concentration  of  the 
pollutant. 


III.  RESULTS  AND  DISCUSSION 
3.1  Spatial  and  temporal  variation  of  pollutants 

Figure  1  represents  the  spatial  evolution  of  the 
concentration  of  the  pollutant  (chromium)  as  a  function  of 
the  depth  x  for  the  different  times  t  =  1  day,  ti  =  5  days  and 
t2  =  30  days,  c/co  is  the  ratio  of  the  simulated  value  of  the 


pollutant  concentration  to  a  depth  x.  The  analysis  of  this 
figure  shows  that  at  a  time  t  the  concentration  of  the 
pollutant  decreases  exponentially  inside  the  filter.  In 
general,  this  concentration  increases  with  the  time. 
Considering  40%  of  the  initial  chromium  concentration  as  a 
guideline,  we  found  that  the  depths  of  the  column  x  should 
be  xi  =  10  cm  for  ti  =  lday,  X2=  20  cm  for  t2=  5  days  and  xj 
=  25,  5  cm  fort3  =  30  days. 

Figure  2  shows  the  time  evolution  of  the  concentration  of 
the  pollutant  at  the  filter  outlet  for  different  depth  values  of 
the  filter:  xi=  20cm,  X2=  25  cm  and  X3=  30  cm  It  is  found 
that  the  concentration  of  the  pollutant  increases  with  time 
and  decreases  as  the  depth  of  the  filter  increases. 


Fig.  1:  Evolution  of  chromium  concentration  in  space  for  x 
=  30  cm 


Fig.  2:  Evolution  of  the  chromium  concentration  overtime 
for  t  =  20  days 
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3.2  Effect  of  the  decay  rate  of  Chromium 


Fig.  3:  the  temporal  evolution  of  the  chromium 
concentration  for  the  values  of  A  =  0.0  and  A  =  0.1  at  t  =  1 0 
min  and  x  =  30  cm. 


Fig.  4:  the  spatial  evolution  of  the  chromium  concentration 
for  the  values  of  A  =  0.0  and  A  =  0.1  at  t  =  10  min  and  x  = 
30  cm. 

Figure  3  represents  the  temporal  evolution  of  the  pollutant 
concentration  in  the  filter  with  different  values  of  the  decay 
rate  X  =  0.0  and  X  =  0.1.  We  find  that  in  the  absence  of  this 
degradation  coefficient  i.e.  when  X  =  0,  the  concentration  of 
the  pollutant  evolves  rapidly  inside  the  filter,  in  the 
presence  of  the  latter,  the  concentration  evolves  slowly. 


ISSN:  2456-1878 

Figure  4  represents  the  spatial  evolution  of  the 
concentration  in  the  filter  with  different  values  of  the  decay 
rate  X  =  0.0,  0.1  and  X  =  2.  The  analysis  of  this  figure 
reveals  an  exponential  decay  of  the  concentration  inside  the 
filter.  The  larger  the  lambda,  the  higher  the  concentration  of 
the  pollutant  decreases  rapidly  inside  the  filter  until 
reaching  the  asymptote  c  =  0.  Therefore,  the  natural  decay 
of  pollutants  plays  an  important  role  in  wastewater 
treatment  since  it  allows  to  avoid  material  waste.lt  is 
necessary  for  the  design  of  industrial  sorption  columns 

3.3  Validation  of  Analytic  and  numerical  approach 


Fig.6:  spatial  variation  of  Concentration 


Fig.  7:  temporal  variation  of  concentration 
Numerical  solutions  of  dispersion  equation  are  generated 
using  the  RK4  algorithm.  For  the  stability  criterion,  we  took 
space  step  and  temporal  step  Ax=  1  and  At  =  1  respectively. 
These  results  were  obtained  fordifferent  values  of  the  filter 
depth  xi  =  20  cm,  X2  =  25  cm  and  X3  =  30cm  and  for  several 
times  ti  =  1  day,  t2  =  5  days  and  t3  =  30  days.  After  analysis, 
it  appears  on  Figure  6that  the  concentration  of  the  pollutant 
decrease  exponentially  with  time  inside  the  water  filter. 
Figure  7  shows  that  the  concentration  of  the  pollutant 
increases  at  the  exit  of  the  filter.  The  variation  of 
theConcentration  of  Pollutants  into  Filter  of  Filter  is  due  to 
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the  concentration  at  the  filter  entrance  and  to  the  1.4.  Spatio-temporal  evolution  of  the  concentration 

characteristics  of  the  material  (ironoxide  coated  sand)  used.  of  the  pollutant 

Finaly,  the  results  obtained  shows  that  these  two  methods 
(analytical  and  numerical)  are  in  good  agreement. 


Figure  8  shows  the  spatio-temporal  evolution  of  the 
concentration  of  the  pollutant  in  the  filter.  The  x-axis  is 
associated  to  the  time  taken  by  the  pollutant  in  the  filter, 
and  the  y-axis  represents  the  depth  of  the  column.  This 
figure  is  subdivided  into  three  phases.  The  initial  phase 
which  starts  from  0  to  30%  of  the  initial  concentration 
characterized  by  the  blue  color,  the  unsaturated  phase  from 
30  to  70%,  and  the  saturated  phase  which  start  from  70  to 
100%  characterized  by  the  red  color.  Figure  8  shows  the 
concentration  of  pollutant  increasewith  time  and  decrease 
with  the  columndepth.  This  results  could  be  used  to  design 
equipment  for  chromium  removal  in  aqueous  solution. 

IV.  CONCLUSION 

The  degradation  coefficient  is  an  important  factor  which 
control  the  design  of  membrane  filter  to  remove  chromium 
VI  in  water.  This  coefficient  must  be  known  to  optimize 
chromium  VI  adsorption  in  aniron  oxide  coated  sand  based 
filter  on  a  technical  scale.  For  adsorption  of  chromium  VI 
on  iron  oxide  coated  sand,  analytical  and  numerical  data 
obtained  from  simulating  a  fixed  bed  adsorber  showed  that 
iron  oxide  coated  sand  was  very  effective  in  removing 
chromium  Wrom  an  aqueous  solution.lt  is  shows  that  the 
degradation  coefficient  of  chromium  X  has  a  significant 
effect  on  the  spatio-temporal  evolution  of  the  concentration 
of  the  pollutant.  This  coefficient  should  be  taken  into 


account  in  models  in  order  to  avoid  to  overestimate  the 
concentration  of  the  pollutant  in  reactive  media.  A  clear 
improvement  is  observed  for  a  material  having  a  high 
adsorption  capacity  and  a  high  first  order  degradation 
coefficient.  Affordable  water  treatment  filters  could  be 
constructed  using  iron  oxide  coated  sand  which  is  an 
effective,  inexpensive  material  for  the  removal  of  water 
pollutants  or  for  the  treatment  of  effluents  from  industries. 
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